Valence bond description of the long-range, nonfrustrated Heisenberg chain 
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The Heisenberg chain with antiferromagnetic, powerlaw exchange has a quantum phase transition 
separating spin liquid and Neel ordered phases at a critical value of the powerlaw exponent a. The 
behaviour of the system can be explained rather simply in terms of a resonating valence bond state in 
which the amplitude for a bond of length r goes as r~ a for a < 1, as r~^ 1+a ^ 2 for 1 < a < 3, and as 
r~ 2 for a > 3. Numerical evaluation of the staggered magnetic moment and Binder cumulant reveals 
a second order transition at a c = 2.18(5), in excellent agreement with quantum Monte Carlo. The 
divergence of the magnetic correlation length is consistent with an exponent v = 2/(3 — a c ) ~ 2.4(2). 



Introduction — Quantum spin-half chains whose inter- 
actions are local and only weakly frustrating^ have a 
quasi-long-range ordered ground state with powerlaw 
spin correlations^ This is different from the situation 
in higher dimensions, where such models exhibit true 
long-range order (LRO)M It is well known that LRO in 
one dimension is proscribed by theorem^ but only when 
the interactions are sufficiently short-ranged. With the 
addition of an antiferromagnetic interaction of arbitrary 
strength and range, the Heisenberg spin chain acquires a 
phase diagram that includes both spin liquid and Neel- 
ordered regions. 

Laflorencie and coworkers^ have proposed a model of 
the form H — Y]^ JijSi ■ Sj with an exchange coupling 
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Here, A and a are positive parameters, is the dis- 



tance between sites i and j, and jij = S(\ 



1) is the 



nearest-neighbour (NN) matrix. The authors of Ref. 6 
have mapped out the X-a phase diagram using quantum 
Monte Carlo. They report the existence of a line of crit- 
ical points — separating the magnetically disordered and 
ordered phases — along which the critical exponents vary 
continuously; the dynamical exponent obeys the inequal- 
ity z < 1. Studies of the A = 1 model have previously 
been carried out using a real-space renormalization group 
method^ and spin wave theory^ 

In this paper, we show that a valence bond (VB) 
description^ of the long-range spin chain provides a uni- 
fied picture of the quantum phase transition and of the 
(seemingly) quite different ground states on cither side 
of it. Moreover, we show that numerical results based 
on a resonating valence bond (RVB) wavefunction are 
quantitatively accurate, even near criticality. 

As we have argued elsewhere ; 10 ' 11 spin models with 
local, nonfrustrating interactions are well described by 
RVB states^ in which bonds of length r appear with a 



probability amplitude (£ + r ) 



-(d+l)/2 



Although this 



result is most accurate for large dimension d, it remains a 
good approximation for spin-5 1 systems in d — 1 provided 
that 2S is odd. The only subtlety is to explain why the 



properties of the 1/r 2 state in d 
in d = 2 are so different. 



1 and the 1/r 3 state 



FIG. 1: A slice in the space of RVB wavefunctions spanned 
by the lattice dimension d and the powerlaw exponent p of 
the long-bond tail. Spin models with local, nonfrustrating 
interactions live along the line p = d + 1. Those with long 
range, nonfrustrating interactions live in the band d < p < 
d + 1. The nearest neighbour Heisenberg model is indicated 
with an open circle (d = 1) and an open square (d = 2). In 
one dimension, the path from p — 2 to p — 1 (from open circle 
to filled circle) passes through a magnetic transition. 



In the VB picture, the liquidness of the linear chain 
and the antifcrromagnctism of the square lattice (for in- 
stance) are consequences of their different lattice geome- 
tries. Overlaps of VB states form a collection of loops, 
and the existence of magnetic LRO is related to the loops 
being macroscopic in size.— In general, the relationship 
between bonds and loops is nontrivial. For the family of 
RVB states with r~ p bond amplitudes, it turns out that 
there is a critical value of the exponent p below which the 
typical loop becomes system-spanning (marking the on- 
set of magnetic LRO). This critical value increases mono- 
tonically with the dimension of the lattice (and diverges 
for d = 3). For the linear chain the value is p c sa 1.6 and 
for the square lattice p c 3. 3 11 ' 14 ; clearly, d+ 1 > 1.6 for 
d = 1 (disorder), whereas d + 1 < 3.3 for d = 2 (order). 

With the introduction of sufficiently long-range inter- 
actions (powerlaw decay with exponent a < 3) , the decay 
exponent of the bond amplitude function can be tuned 
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r' + r" - 1 



FIG. 2: A valence bond is a singlet constructed from two spins 
in opposite sublattices and characterized by the directed dis- 
tance between them. In a valence bond state, all the spins 
are paired to form valence bonds in one of the {L/2)\ possible 
configurations. The reconfiguration induced by the nearest- 
neighbour Heisenberg interaction is sketched above. The dot- 
ted (red) line shows where the interaction is applied. Bonds 
r',r" are mapped to 1,7-' + r" — 1. 



continuously in the range 1 < p < 2. See Fig. [T] For the 
linear spin chain (in a model related to Eq. ([1]) with A = 
1), the critical p c is achieved at a c = 2.18(5). This com- 
pares favourably to the critical value a^ mc = 2.225(25) 
determined by quantum Monte Carlo in Ref. [6] and repre- 
sents a significant improvement over the values predicted 
by spin wave theory, a™ = 2. 46, 6 and by a numerical 
renormalization group method, a r c s = 1.85.— 

RVB analysis — The singlet ground state of an even 
number of S = ^ spins can be expressed in an overcom- 
plete basis of bipartite VB states*^ A simple RVB wave- 
function can be constructed by assigning to each bond 
of odd length r an amplitude h(r) and taking as the to- 
tal weight for each VB configuration the product of the 
amplitudes of the individual bonds: 



i*>=e n*( 
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(2) 



The sum in Eq. ||3J| is over all possible pairings of spins 
in opposite sublattices, and the product is over all bond 
lengths r appearing in the VB configuration. 

The proper choice of h(r) is determined by the model 
at hand. We observe that the application of a bipartite 
Heisenberg interaction to a VB state results in the re- 
configuation of bonds depicted in Fig. [2l When acting 
between sites i and j (in opposite sublattices), separated 



by a distance a 



j , the interaction transforms bonds of 



length r' and r" into bonds of length a and r = r'+r" — a. 
The steady-state solution-^ of this reconfiguration pro- 
cess is 



h(r) 
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where J q is the Fourier transform of the interaction nor- 
malized to J q= o = 1- In a model with NN interactions 
only, J q = cos g, and the long distance behaviour is 
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FIG. 3: These log-log plots illustrate the five different scaling 
behaviours described in Eq. (|10[) . Various renormalizations of 
the bond amplitude h are plotted with respect to the reduced 
coordinate p — 2r/L. All the vertical axes (labels suppressed) 
should be understood to start at 10° and increase by one 
decade per tic (10 1 , 10 2 , etc.). 



Numerical evaluation of Eq. ([2|) can be performed 
stochastically on lattices of finite size, as described in 
Ref. M 

We find that the RVB state for the Heisenberg 
chain (characterized by Eq. ([3]) with J q — cosg) has a 
magnetically disordered ground state and an L — > oo ex- 
trapolated energy E = -0.4360(1), within 1.7% of the 
exact result E = log 2 - 1/4 = -0.44315. The dis- 
crepancy is somewhat large, but not unreasonably so 
given that our wavefunction was not variationally deter- 
mined. Moreover, the RVB state performs worst in the 
disordered phase; the agreement is increasingly good the 
deeper we go into the magnetic region. 

We now consider the long-range exchange integral 
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which is equivalent to the A = 1 case in Eq. ((T]), ex- 
cept that we have removed the interactions between spins 
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FIG. 4: The long-bond behaviour of the amplitude func- 
tion is powerlaw irrespective of the range of interaction, but 
the decay exponent changes as a function of a. We identify 
three regimes. When a > 3, the bond amplitudes share the 
h(r) ~ r~ 2 and uj q ~ q behaviour of the nearest-neighbour 
Heisenberg model, differing only in the amplitudes of the 
short bonds and the value of the spinon velocity. In the 
intermediate regime 1 < a < 3, the decay exponent varies 
continuously, and the spin waves exhibit sublinear dispersion, 

LOq ~ q 



(« !)/ 2 Below a critical (A-dependent) a c , the spins 



develop long-range antiferromagnetic order. When a < 1, the 
interaction is superextensive and the system is close to being 
classically Neel ordered. 



in the same sublattice and compensated with coupling 
strength of opposite sign at neighbouring sites. This 
change has no real significance but it does simplify our 
analysis since nonbipartite interactions would require in- 
troducing a second update rule (different from the one 
shown in Fig. (2). 

In this case, the J q to appear in Eq. ([3]) is 
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Here, n runs over all odd integers from 1 to L/2 — 1, and, 
to leading order in L, 
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When a > 3, £ Q is 0(1), and the long wavelength be- 



haviour of h„ 



h{r) 



0(^q 3 ) is integrable; hence, 
2^ 



This is the same form we found for J q = cos q, except 
that £ Q , which controls the length scale at which the r~ 2 
tail is cut off, is now a dependent. In the limit a — > oo, 
£ a — » 1, and we recover Eq. (g]). 

As a — > 3 + , ^ a diverges, and h{r) changes functional 
form to r _2 (logr) 1 / 2 at a = 3. There is an additional 
crossover at a = 1, where h(r) goes as r _1 (logr) -1 / 2 . 
Otherwise, the bond amplitude has a continuously vari- 
able decay exponent: r~ a for cv < 1, and r~( 1+a )/ 2 for 
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FIG. 5: (M 2 ) = (l/£ 2 )E«(-l) i+3 '<Si ■ S,), the staggered 
structure factor, is computed for lattices up L — 8192 in size. 
The main plot shows the magnitude of the staggered magnetic 
moment, M = (M 2 ) 1 ^ 2 , as a function of a. The value of M in 
the thermodynamic limit is marked with a solid (green) line. 
The inset shows the finite-size scaling of (M 2 ) in the vicinity 
of the magnetic quantum critical point. The extrapolations 
to 1/L = are draw with solid (green) or dashed (orange) 
lines depending on whether the intersection with the y-axis is 
positive. We estimate that the critical point lies in the range 
2.1 < a c < 2.2. 



1 < a < 3. These results can be demonstrated by noting 
that h(r) obeys the relation 
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with a scaling function 



H a (p) 
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as verified by data collapse in Fig. [3l 

From what we know about the RVB ground state, it 
is possible to make a reasonable guess for the excitation 
spectrum u> q . Wherever the system is Neel ordered, the 
low-lying excitations are magnons. These spin-1 exci- 
tations can be created by promoting in succession each 
singlet bond in the RVB state to a triplet (along with the 
appropriate ^-dependent phase factor): at the mean field 
level, one finds that h q ~ e~" 9 in the q — > limit.— With 
this identification, u> q ~ q when a > 3, and uj q ~ g( Q_1 )/ 2 
when 1 < a < 3. 

Numerical evaluation of the RVB wavefunction over 
a range of a reveals the phase diagram summarized in 
Fig. 2J We find that the staggered moment, shown in 
Fig. is nearly saturated at M = 0.5 throughout the 
a < 1 semi-classical region. M decreases monotonically 
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FIG. 6: (Left) Q = 1- (Mj >/3{Aff } 2 , the Binder cumulant, is 
plotted as a function of a for lattices of increasing size. (The 
curves have increasing slope, approaching a step function for 
L — > oo). The data points are shaded as in Figs. [3] and [S] The 
lines, connecting points of equal L, emphasize the invariance 
of Q at q = a c . (Right) The same data are replotted with an 
a c shift and L 1 ^ rescaling of the horizontal axis. Good data 
collapse is acheived for a c — 2.18 and v — 2/(3 — a c ) — 2.44. 
The line (1/5) [1 — tanh(a;/5)] is drawn as a guide to the eye. 

across the intermediate regime 1 < a < 2 and vanishes 
continuously at some a — a c . Quantum disorder for 
a > 3 is guaranteed by theorem^ 

Under the assumption of Lu q ~ gt"" 1 )/ 2 S pi n excita- 
tions, the staggered structure factor exhibits 1/L^ 3_Q ^ 2 
scaling.S. Accordingly, the inset in Fig. [5] shows it van- 
ishing somewhere in the range 2.1 < a c < 2.2. A 
more sophisticated finite-size scaling analysis (see Fig.[B]), 
based on data collapse of the Binder cumulant^ indi- 
cates a c — 2.18(5). The correlation length exponent ex- 
tracted from the fit, v — 2.4(2), is in agreement with the 
large-A prediction of the corresponding A-component 
vector theory £ namely, v = l/(a c — 1) for 1 < a c < 5/3, 



and v — 2/(3 — a c ) for 5/3 < a c < 3. 

The transition is in the quantum percolation classic 
As the range of interaction is ramped up, the bonds in 
the RVB state grow longer, and the overlap loops in- 
crease in size. Across the magnetic transition, the scal- 
ing dimension D of the average loop size (~ L D ) changes 
discontinuously from D = to D = 1. At criticality, the 
loops have fractal dimension < Df < 1, and the dy- 
namical exponent is fixed by z = Df. This explains the 
observation in Ref.0that z < 1. (Analogous behaviour- 1 
is found in two dimensions for radially symmetric bond 
amplitude functions.) 

Conclusions — The NN Heisenberg chain has a quan- 
tum disordered ground state, but the addition of long- 
range, antiferromagnetic interactions can drive the for- 
mation of magnetic LRO. The ground state wavefunction 
evolves smoothly across the phase boundary. Its struc- 
ture on both sides of the transition is essentially that of 
a factorizable RVB wavefunction whose bond amplitudes 
decay as a powerlaw in the bond length. Only the value 
of the decay exponent changes across the transition. 

The onset of Neel order can be understood as a quan- 
tum percolation transition in which the VB loops become 
system-spanning. The dynamical exponent < z < 1 is 
equal to the fractal dimension of the VB loops that are 
formed at criticality. The exact point at which the tran- 
sition occurs and the value of the critical exponents there 
depend sensitively on the bond distribution. 

The RVB picture is often invoked in a loose, heuristic 
way. Here, we have presented a concrete RVB wavefunc- 
tion that proves to be a practical computational tool. 
Many aspects of the unbiased quantum Monte Carlo re- 
sults can be reproduced by the RVB wavefunction with 
dramatically less computational effort. 
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